function out = int_jt(t, r)
%INT_JT Summary of this function goes here
%   Detailed explanation goes here

global q;
global A;
global h_jt;

tmp_phi = phi(t, r);
diff_phi_t = (tmp_phi - phi(t - h_jt, r)) ./ h_jt; 

out = imag(conj(tmp_phi) .* (diff_phi_t - 1i .* q .* A(r)) .* tmp_phi);

end

